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Abstract 

QCD sum rules are based on the Operator Product Expansion of current correlators, 
and on QCD-hadron duality. An extension of this program to finite temperature is 
discussed. This allows for a study of deconfinement and chiral-symmetry restora- 
tion. In addition, it is possible to relate certain hadronic matrix elements to ex- 
pectation values of quark and gluon field operators by using thermal Finite Energy 
Sum Rules. In this way one can determine the temperature behaviour of hadron 
masses and couplings, as well as form factors. An attempt is made to clarify some 
misconceptions in the existing literature on QCD sum rules at finite temperature. 
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1 Outline 



This report is an expanded version of a talk given at CAM-94, the joint meeting of 
the Canadian Association of Physicists, the American Physical Society, and the Mexican 
Physical Society. I discuss here an extension of the QCD sum rule program to finite 
temperature. At T = 0, QCD sum rules offer a very successful quantum field theory 
framework to extract information on hadronic physics from QCD analytically. This great 
success, unfortunately, does not translate immediately into the finite temperature domain, 
where there are some serious unresolved problems with Laplace transform sum rules. In 
spite of this, some reasonable progress has been made through lowest moment Finite 
Energy Sum Rules (FESR). As with any new field, some degree of confusion is to be 
expected as our ideas take shape, and we gain a better understanding of the subject. It is 
important, though, not to persist on wrong notions, nor advocate results which have been 
shown to be in contradiction with more fundamental facts. 1 shall attempt to clarify here 
some misconceptions present in the existing literature on this subject. The purpose is not 
to antagonize, but rather to point out the problems and promote discussion that might 
help to solve them. I may group these misconceptions into two categories: (a) conceptual, 
and (b) specific results of applications. Among the first is the statement (without proof) 
that the notion of QCD-hadron duaUty, one of the two pillars of QCD sum rules, should 
abruptly disappear as soon as the temperature is turned on (by no matter what small 
amount). If correct, this proposition would invalidate a smooth extension of the QCD 
sum rule program to T 7^ 0. I shall argue against this scenario. Second, and in connection 
with results of applications, there exist quite a few determinations of hadron masses based 
on Laplace transforms which ignore the underlying problems with these sum rules. Chief 
among these problems is the fact that they are inconsistent with the well known (and well 
established) temperature behaviour of the quark condensate and the gluon condensate. 
Not surprisingly, the predictions from this approach are in serious contradiction with the 
T-dependence of hadron masses obtained in other independent frameworks. 

The outline of this report is as follows. In Section 2 I introduce briefly the key ideas 

behind QCD sum rules at T = 0, and argue for a smooth extension of this program to 
finite temperature. 1 discuss supportive evidence for the validity of both the Operator 
Product Expansion, and the notion of QCD-hadron duality at T 7^ 0. Section 3 deals with 
chiral-symmetry restoration, and in Section 4 I review proposals for phenomenological 
order parameters to characterize deconfinement. In Section 5 I show how these two phase 
transitions have been related using a lowest moment FESR-QCD Sum Rule. Section 
6 summarizes recent results for the temperature behaviour of the electromagnetic form 
factor of the pion in the space-like region, once again through a lowest moment FESR. 
Finally, in Section 7 I concentrate on the problems with Laplace transform QCD sum 
rules, as well as with higher moment FESR, and make critical comments on the existing 
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literature. Except for a major part of Section 7, the ideas and results presented here have 
already been published in the literature, and the reader may trace them through the list 
of references. 



2 QCD sum rules 



QCD sum rules at T = |I| are based on the Operator Product Expansion (OPE) 
of current correlators at short distances, suitably extended to include non-perturbative 
effects. The latter are parametrized in terms of a set of vacuum expectation values 
of the quark and gluon fields entering the basic QCD Lagrangian. Contact with the 
hadronic world of large distances is achieved by invoking the notion of QCD-hadron 
duality. This leads to relationships between fundamental QCD parameters {Kqcdi quark 
masses, vacuum condensates, etc.) and low energy parameters (hadron masses, widths, 
couplings, form factors, etc.). The values of the vacuum condensates in the OPE cannot be 
calculated analytically from first principles, as this would be tantamount to solving QCD 
exactly. Instead, they are extracted from certain channels where experimental information 
is available, e.g. e"'"e~ annihilation, and r decays p. It is also possible, in principle, to 
estimate them numerically from lattice QCD. To be more specific, let us consider the 
two-point function 

Y[{q)=i j (fx exp(igx) < 0|T(J(x), J"^(0))|0 > , (2.1) 

where J{x) is a local current built from the quark and/or gluon fields entering the QCD 
Lagrangian, and having specific quantum numbers. In the sequel I concentrate on light 
quark flavours. The OPE of n(g) is formally written as 

U{q) = Ci<I> <Or>, (2.2) 

r 

where the Wilson coefficients Cr{q) depend on the Lorentz indices and quantum numbers 
of the external current J{x), and also of the local gauge-invariant operators Or built from 
the quark and gluon fields of QCD. The unit operator I in Eq.(2.2) represents the purely 
perturbative piece. The OPE is assumed valid, even in the presence of non-perturbative 
effects, for < (spacelike), and | |^ ^qcd- principle, all Wilson coefficients are 
calculable in perturbative QCD to any desired order in the strong coupling constant. In 
the sequel we shall work at leading (one loop) order for simplicity. The non-perturbative 
effects are then buried in the vacuum condensates. Since these have dimensions, the 
associated Wilson coefficients fall off as inverse powers of = —q^. For instance, if the 
current J(x) in Eq.(2.1) is identified with the axial- vector current A^{x) =: u{x)'~^^'~^^d{x) :, 
then with 

n^.(g) = -Qt^uMQ^) + q>.qMq^) , (2.3) 
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one easily finds 

47r no(g) = - In — + + — + ■■■, (2.4) 

where /i is a renormalization scale, and e.g. the leading vacuum condensate is given by 

C4 < O4 >= - < as G'^ > -Sn'^Thg < qq > , (2.5) 
3 

with rhq = {rriu + md)/2, and < qq >=< uu >~< dd >. In Eq.(2.4) a term proportional 
to m? JQ"^ has been neglected. The function Iio{q), Eq. (2.4), satisfies a dispersion relation 
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defined in this case up to one subtraction constant, which can be disposed of by e.g. taking 
the first derivative with respect to in Eq.(2.6). The notion of QCD-hadron duality 
is implemented by calculating the left hand side of Eq.(2.6) in QCD through the OPE, 
and parametrizing the spectral function entering the right hand side in terms of hadronic 
resonances, followed by a hadronic continuum modelled by perturbative QCD. In this 
fashion one relates fundamental QCD parameters, such as quark masses, renormalization 
scales, vacuum condensates, etc., to hadronic parameters such as particle masses, widths, 
couplings, etc.. The convergence of the Hilbert transform, Eq.(2.6), may be improved by 
considering other integral kernels. This leads to other versions of QCD sum rules, such as 
the Laplace transform. Finite Energy Sum Rules (FESR), etc.. For instance, the Laplace 
transform QCD sum rule for the axial-axial correlator is 

/•oo 1 

C no(Q2) = no(M2) = / ds exp(-s/M2) - Imli^is) , (2.7) 

Jo TV 

where is the Laplace parameter which plays the role of as the short distance 
expansion parameter (for light quark correlators). Performing the Laplace transform of 
Eq.(2.4) one finds 

i -M-s/M') - /.nn,(.) = (1 + ^jpi- + ^ + ■■■). (2.8) 

One should notice from Eq.(2.8) that in Laplace transform QCD sum rules all condensates 
are involved. Their numerical importance, though, is suppressed by inverse powers of M^, 
as well as by factorial coefficients. 

Either by using Cauchy theorem, or by expanding Eq.(2.7) in M^, the Laplace transform 
QCD sum rule Eq.(2.7) is equivalent to the infinite number of FESR 

i-f-' C2N < O2N >= 47r2 / ds s""-' - ImU,is)\nES - % , (2.9) 

JO TT iV 
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where N=l,2,..., and I'mIlQ{s)\fiES stands for the purely resonant contribution to the 
hadronic spectral function (i.e. without the continuum). For = 1, C2 < O2 > is nothing 
but a short-hand notation for the perturbative quark mass insertion (C2 < O2 >oc m'^). 
A practical advantage of these FESR is that now the vacuum condensates of different 
dimensionality are effectively decoupled. This becomes particularly important if the ob- 
jective is to extract the values of these condensates from a knowledge of the hadronic 
spectral function. In fact, this is the correct procedure followed in [0] to determine the 
condensates from data on e~^e~ annihilation and r-decays. The decoupling of the different 
Ctv < Oat > is also important to study the self-consistency of the Laplace transform QCD 
sum rules at T 7^ 0, as will be discussed in Section 7. 

An extension of this QCD sum rule program to finite temperature was proposed some 
time ago in 0]. This proposal entails the assumptions that (a) the OPE continues to 
be valid, except that now the vacuum condensates will develop an {a-priori) unknown 
temperature dependence, and (b) the notion of QCD-hadron duality also remains valid. 
I shall discuss below some evidence in support of these two assumptions. Notice that in 
analogy with the situation at T = 0, the thermal behaviour of the vacuum condensates is 
not calculable analytically from first principles. Some model or approximation must be 
invoked, e.g. the dilute pion gas approximation, lattice QCD, etc.. The quark, the gluon, 
and the four-quark condensates at T 7^ have thus been estimated in the literature [^]- 
0. At finite temperature, the basic object to be considered is the retarded (advanced) 
two-point function after appropriate Gibbs averaging 

n (g, T) = z J d^x exp (iqx) 9{xo) « [J(x), J\0)] » , (2.10) 

where 

«A-B »= J2 exp{-Er,/T) {n\A ■ B\n) /Tr{exp{-H/T)) , (2.11) 

n 

and I n > is a complete set of eigenstates of the (QCD) Hamiltonian. The OPE of n(g, T) 
is now written as 

n (g, T) = Ci « I » + J2 Cr{q) « Or » , (2.12) 

r 

It must be stressed that the states | n > entering Eq.(2.11) can be any complete set of 
states, e.g. hadronic states, quark-gluon basis, etc.. The hadronic (mostly pion) basis 
has been advocated in [Q, while the quark-gluon basis was first used in |^. These two 
approaches are quite complementary, rather than in conflict, as the information they 
provide is somewhat different. The pion basis is well suited to determine the temperature 
dependence of vacuum condensates at low T. It does not make use of QCD-hadron 
duality, and thus has little relationship to the QCD sum rule program. On the other 
hand, use of the quark-gluon basis allows for a smooth extension of that program to 
finite temperature. As it continues to rely on both the OPE and QCD-hadron duality, 
this approach provides information on thermal Green functions provided the temperature 
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dependence of the condensates is known. Since the latter can be obtained e.g. from using 
the pion basis in Eq.(2.11), the two choices of the complete set | n > complement each 
other. However, it must be kept in mind that the choice of the pion basis, being a form 
of the virial expansion, is restricted to low temperatures. This is not necessarily the case 
for the quark-gluon basis approach. In fact, given an expression for the condensates in 
some framework, accurate enough for all T up to the critical temperature, the QCD sum 
rules will provide the T-dependence of hadronic matrix elements in the same temperature 
range. 



The validity of the OPE (at any temperature) beyond perturbation theory cannot be 
proven from first principles, since one does not know how to solve QCD exactly. However, 
at T = one can solve exactly other field theories which bear some resemblance to 
QCD, thus providing evidence in support of this assumption. For instance, a study of 
the (exactly solvable) 0{N) sigma model, in the large N limit, and of the Schwinger 
model, both in two dimensions, shows that the short distance approximation to exact 
Green functions agrees with the result from the OPE [Q. An extension of this analysis 



to finite temperature [|T0| shows the same agreement, and thus supports the assumption 



of the validity of the OPE in this regime. I briefly summarize the results of |]T0[- Let 
me consider first the 0{N) sigma model in 1+1 dimensions which is characterized by the 
Lagrangian 

C = ^[d,a^{x)][d,a'^{x)], (2.13) 

where a = 1,...N, and a°'a°' = N/f, with f being the coupling constant. In the large N 
limit this model can be solved exactly (for details see 0), it is known to be asymptotically 
free, and in spite of the absence of mass parameters in Eq.(2.13), it exhibits dynamical 
mass generation. In addition, in this model there are vacuum condensates, e.g. to leading 
order in 1/N : (0|a|0) = \/N m? , whereas all other condensates factorize, viz. (0|a'^|0) = 
{\/N rn?Y ■ The a field is: a = /(9^(T")^/A/iV, and we are interested in the Green 
function associated with the propagation of quanta of this a field. We have calculated 
this Green function at finite temperature [1^. Its imaginary part can be integrated 
analytically in closed form and is 

Im T{uj, q = 0, T) = [1 + Sn^ (cu/2T)] 
2a; 



1 

+ 2 



<< a » 6 << >> 



where the first term above corresponds to the perturbative contribution, the second to 
the non-perturbative, and ub is the thermal Bose factor. Equation (2.14) is valid in the 
time-like region; the space-like region counterpart vanishes in 2 dimensions. Since the 
model is exactly solvable, the thermal behaviour of the vacuum condensates can also be 
calculated, viz. 

« a »=< a > [l + 3nBiuj/2T)] (2.15) 
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In this case the vacuum condensates contribute to the imaginary part, and as Eq.(2.14) 
shows, the thermal dependence of the perturbative piece cannot be absorbed into the 
condensates. Hence, no confusion should arise between perturbative and non-perturbative 
contributions. 

Next, I consider the Schwinger model in 1+1 dimensions, with the Lagrangian 

C = -^F^,F^, + ^i-i^V^ij (2.16) 

where = id^ + eA^. At T=0 this model has been solved exactly, and in the framework 
of the OPE [P]. The short distance expansion of the exact solution coincides with that 
from the OPE. Here, we are interested in the two-point functions 

H++(x) = (0|r{j+(x)j+(0)}|0) (2.17) 

H+_(x) = (0|r{j+(x)j-(0)}|0) (2.18) 

where the scalar currents are: j"*" = i'R ipi , j~ = V'l i^R , with 4'l,r = (1 =t 75)V^/2. 
The function H++((5) vanishes identically in perturbation theory, and the leading non- 
perturbative contribution involves a four-fermion vacuum condensate. We have calculated 



the thermal behaviour of these current correlators [10| and obtain, e.g. for their imaginary 



parts in the time-like region (again, there is no space-like contribution in 2 dimensions) 

ImH++(tu,q = 0,T) =0 (2.19) 

Im n+_(cj, q = 0, T) = ^ [1 - 2n^(cj/2T)] (2.20) 

Hence, the choice of the fermion basis in the Gibbs average of current correlators does not 
imply confusing these fermions with condensates, as argued in p]. As Eqs.(2.19)-(2.20) 
indicate, (perturbative) fermion loop terms and (non-perturbative) vacuum condensates 
develop their own temperature dependence, which in this particular example happen to 
be different. 



Concerning the notion of QCD-hadron duality, it has been suggested recently that it 
is not applicable at finite temperature. If correct, this would require a singular dynamical 
mechanism of a discontinuous nature in order to invalidate the inter-relationship between 
QCD and hadronic parameters effected by duality. No such mechanism has been proposed 
in p|. That this inter-relationship would abruptly disappear by raising the temperature 
from T = to some arbitrary small value, say a nano-Kelvin, seems quite unlikely, es- 
pecially in the absence of a concrete mechanism to achieve it. According to the QCD 
sum rule philosophy, at T = one calculates the theoretical left hand side of Eq.(2.6) 
through the OPE Eq.(2.2), i.e. one uses quark-gluon degrees of freedom, and duality 
relates this QCD part to a weighted average of the hadronic spectral function. The latter 
arises from using hadronic degrees of freedom. At very low temperatures the hadronic 
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spectrum is expected to change very little, and the external current will still convert into 
quark- antiquark pairs. The temperature dependence of the quark and gluon condensates 
is known, and at very low T they also hardly change. Hence, it is only reasonable to 
assume that nothing drastic will happen to duality. There is a sort of temperature inertia 
affecting both QCD and hadronic physics at very low T. At finite temperature, though, 
there are some new effects coming into play, e.g. there are contributions to the QCD 
and hadronic spectral functions in the space-like region (as opposed to only the time-like 
region at T = 0), and the heat bath can support condensates with non-trivial quantum 
numbers. However, these additional contributions vanish smoothly as T approaches zero, 
i.e. they do not introduce any discontinuous behaviour that would abruptly invalidate 
the notion of QCD-hadron duality. At moderate temperatures the hadronic spectrum 
is expected to suffer some rearrangement, in pace with changes in the condensates and 
the increasing importance of the new analytic structure in the complex energy plane. 
By retaining the notion of QCD-hadron duality one is able to relate quantitatively the 
temperature dependence of hadronic parameters with that of QCD parameters, as will be 
discussed in Sections 5-6. 



3 Chiral-Symmetry and its restoration 



1 begin by discussing some of the symmetries of the QCD Lagrangian 

Cqcd = D,-M),j;-^ F,, F^'' , (3.1) 

where 
with 

(a = l,---8) (3.3) 

and M the quark mass matrix. For the purposes of this talk 1 consider only two quark 
flavours (up and down). Among the various symmetries of Cqcv one finds a (global) 
SU (2)y and an SU (2)^^ ® SU (2) ^ symmetry. These Lagrangian symmetries are explicitly 
broken by the quark masses. In fact, the vector (I-spin) Noether current — ■07^ rV, 
and the axial- vector current ^^ = ■0 7^ 75 have divergences 

d^'V^ = i{md - mu)du , (3.4) 

d^A^j^ = i{md + mujd'j^u . (3.5) 
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In the limit rriu = md = 0,SU{2)v, and SU{2)fi (g) SU(2)l become exact Lagrangian 
symmetries. However, this is not what is usually meant by chiral-symmetry restoration, 
which refers to the symmetry of the vacuum. Given a Lagrangian symmetry one must 
investigate how it is realized in the states, starting with the vacuum. According to whether 
the Noether charges = J d^x Jq{x, t) annihilate the vacuum or not, one has a Wigner- 
Weyl or a Nambu-Goldstone realization of the Lagrangian symmetry. In the former 
case particles are classified according to the irreducible representations of the symmetry 
group, as in e.g. SU{2)y. The Nambu-Goldstone realization (spontaneous symmetry 
breaking) corresponds to a hidden symmetry, as the vacuum does not share the symmetry 
of the Lagrangian. This is the case for SU (2)^ ® SU (2) i since e.g. there are no (quasi) 
degenerate parity doublets in the particle spectrum. No particle classification is possible 
in this phase as the vacuum is contaminated by an arbitrary number of massless (Nambu- 
Goldstone) bosons carrying non- trivial quantum numbers. In the case of SU {2)ji^SU{2) l 
the three emerging Nambu-Goldstone bosons are readily identified with the pion (vr^, 7r°), 
which decays to the (hadronic) vacuum through the axial-vector current, i.e. 

(0|A;(0)|7r^■(p))=^^p,5^^ (3.6) 

with = 93.2 MeV. In the limit m„ = = the axial-vector current is strictly 
conserved and, hence, /i^ = 0. In the Nambu-Goldstone phase 

oc (m„ + nid) , (3.7) 

oc (0|uu + J(i|0) ^ , (3.8) 

with the proportionality constants such that 2/|/i^ = (m„ -|- md){Q\uu + dd\Q). A phase 
transition from a Nambu-Goldstone to a Wigner-Weyl mode (chiral-symmetry restora- 
tion) is characterized by the vanishing of the order parameter /7r,or alternatively (0|gg|0). 
Clearly, this can only happen at finite temperature, and if the phase transition does take 
place this should happen regardless of whether the quark masses are zero or not. The 
numerical value of the critical temperature, though, is expected to depend on this fact. 



The temperature behaviour of fniT) at low T {T <^ fi^,) has been investigated in chiral 
perturbation theory with the result M 



U{T) = MO) 



1 - 



^2 



+ 0{T^) 



(3.9) 



in the chiral limit for three flavours. Chiral-symmetry breaking corrections to Eq.(3.9) 
have been also calculated [^], together with corrections due to massive states 0. The 
low temperature expansion of {qq)T has been carried out up to order 0{T^) in ||^. As 
the authors of have pointed out, this low temperature expansion should not be 

extrapolated to the critical temperature, as it is not valid there. For instance, if the phase 
transition is of second order then the critical exponent should be 



\, rather than 2 as one 
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would naively obtain from Eq.(3.9). If one is interested in the behaviour of fn{T) near 
T = Tc, other methods should be used, e.g. the composite operator formalism of |0, 
which reproduces Eq.(3.9) at low T, but gives instead fn(T) oc (1 — T/Tc)^ as T — > Tc. 
This feature will become particularly important in Section 5. 



4 Quark Deconfinement 



At zero temperature the shape of a typical hadronic spectral function consists of some 
delta functions plus resonances with increasing widths, followed by a smooth continuum 
starting at some threshold energy Eq. It is known from fits to actual data that for E>Eq 
the hadronic spectral functions are well approximated by perturbative QCD; Eq is thus 
called the asymptotic freedom (A.F.) threshold. This picture is well supported by all 
exisiting experimental data. Under the assumption that quark deconfinement does take 
place at some critical temperature T^, one would expect that by increasing T from T = 
the resonance peaks in the spectral function should become broader. At T = the reso- 
nance widths would then become infinite, signalling quark deconfinement. This resonance 
melting with increasing temperature would be accompanied by a shift of the asymptotic 
freedom threshold Sq towards threshold. In this picture the resonance width, and/or the 
asymptotic freedom threshold provide a suitable order parameter. 



In the complex s-plane, bound states (e.g. the pion) correspond to poles of the S-matrix 
lying on the real axis, and resonances correspond to poles located in the second Rie- 
mann sheet, their distance to the real axis being measured by the width Tr. At finite 
temperature the form of the Green function will be 

GiE, T) oc -4 ■■ . (4.1) 

In the picture described above, an increase in the temperature will shift all poles farther 
away from the real axis as Tfi{T) increases. While the mass Mr may depend on T, 
it is not a relevant order parameter. With the second Riemann sheet poles infinitely 
far from the real axis, and the A.F. threshold correspondingly close to the origin, the 
spectral function is then described entirely by the QCD continuum extrapolated down to 
the kinematic threshold. One may interpret such a phase as a deconfined phase where all 
hadrons originally contributing to the spectral function have melted into quarks or quark- 
antiquark pairs. As an example, let us consider the following ansatz for the rho-meson 
width 
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where Td is the critical temperature for deconfinement, and a a positive number. The 
behaviour of the hadronic spectral function in the vector channel, with Mp{T) = Mp{T = 
0), Tp as in Eq.(4.2) with a = 0.5, is illustrated in Fig. 1 for two values of the temperature: 
T = (sohd curve), and T = Td/2 (dashed curve). 

In the case of stable hadrons, e.g. the pion or the nucleon, for which r(T = 0) = 0, we 
would expect an imaginary part in their propagators to develop for T 7^ 0. In this case, 
these hadronic thermal widths should be interpreted as damping coefficients of wave pack- 
ets propagating through a dispersive medium (heat bath). Their growth with increasing 
temperature implies that, with hadrons corresponding to excitations in this medium, these 
excitations become less and less important. Combined with the notion of QCD-hadron 
duality, this may be interpreted as the melting of hadrons into their constituents, thus 
signalling a deconfinement phase transition. 



I wish to stress that the relevant order parameter for deconfinement should be the width 
and not the mass. Claims have been made occasionally in the literature that the mass 
should vanish at the critical temperature. While this may happen in some cases, it is 
not a necessary condition for deconfinement. In fact, let us consider a stable hadron, 
i.e. one with zero width at zero temperature, such as the pion or the nucleon. If with 
increasing temperature the mass goes to zero, and nothing happens to the width, then 
at the critical temperature one would still see a peak in the hadronic spectral function 
at zero energy. The hadron has not disappeared from the spectrum, and hence has not 
melted! The only way a particle can melt is by having a temperature dependent width 
such that F/j — >• cxo as T — >• Tc. What happens to the mass (defined as the position of 
the pole on the real axis) is irrelevant to this argument. Once the resonance becomes 
infinitely broad, the spectral function becomes smooth and should be well approximated 
by the quark degrees of freedom, i.e. by perturbative QCD. Since this proposal was first 
made in independent supportive theoretical evidence has become available. In it 
has been shown in the framework of the virial expansion that: (a) at low temperatures 
(T < 50 MeV) and in the chiral limit (yU^ = 0) the mass of the nucleon Mj\f{T) ~ M7v(0), 
and its width T^lT) ~ FAr(O) = 0; (b) At higher temperatures, and now away from the 
chiral limit, both M^iT) and fJ^niT) increase slightly (by 4% and 1%, respectively, up 
to T<160 MeV), while F7v(T) and F^(T) increase substantially, e.g. at T = 160 MeV 
the ratio of width to mass (imaginary to real part of the propagator) is about 43% for 
the pion and 20% for the nucleon. Since the pion and nucleon widths are strictly zero at 
T = 0, this is quite a dramatic effect. Additional evidence for the approximate constancy 
of Mat and /i,r follows from the sigma model [|l5l, which also gives increases over Mjv(O) 
and fj,n{0) at the level of a few percent. Further support comes from a recent calculation 
of ^Tr{T) in the composite operator formalism [0, showing approximate constancy of the 
pion mass over a wide range of temperatures, with a tendency to increase near the critical 
temperature. QCD sum rules give a similar T-dependence of the nucleon mass The 
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temperature behaviour of the pion and nucleon widths have also been calculated recently 



in the framework of the sigma model ||T8[. The results are in qualitative agreement with 
, and are shown in Figs. 2-3 for three different values of the a-meson mass: M„ = 400 
MeV (a), 600 MeV (b), 800 MeV (c). Concerning the rho-meson mass, it has been shown 
19| that to first order in the virial expansion, unitarity of the vr — vr scattering amplitude 



requires that if Tp{T) increases with T, then Mp(T) must also increase with T. I shall 
come back to this point in Section 7. 



5 Relationship between deconfinement and chiral- symmetry 
restoration 



Let us consider the retarded two-point function Eq.(2.10) involving the axial- vector cur- 
rent Ap_{x) =: 'u(a;)7^75(i(x) :. As discussed in Section 2, the QCD sum rule program 
consists in calculating the two-point function in the deep euclidean region in perturba- 
tive QCD, adding non-perturbative effects parametrized in terms of vacuum expectation 
values of the quark and gluon fields appearing in the QCD Lagrangian, and relating the 
result to the hadronic spectral function by means of a dispersion relation. 

In order to find the perturbative QCD behaviour of Eq.(2.10) we assume a dilute quark 
gas at temperature T with zero chemical potential. The virtual quanta associated to the 
local current A^{x) will convert into qq pairs for q^ = u"^ — > ^im?^, while for space-like 
momenta (a;^ — q^) < there is an additional cut in the complex tu-plane centered around 
a; = 1^. These two distinct processes contribute to the spectral function as follows 



-Imn+ (u;, q) = E / LIPS{uj, q, Ei, pi^a, P2) 

X {Q\A^\qq){qq\A,\Q) x [1 - rip^E^) - np^E^)] , (5.1) 

for > 0, and 



q q 

X {q\A^\q){q\A,\q) X [n^(^i) - 72^(^2)] (5.2) 
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for < 0. In the above equations nF{E) is the Fermi-Dirac distribution function, and 
the phase space is 

TTPQf 7? TP \ ^^Pi d^P2 

LIPS{., q, E„ p„ E,, p,) = ^^-^^ ^^-^ 

X 5(cu - El - E2)5^'Hq - Pi - P2) . (5.3) 
It must be stressed that the above two contributions to the spectral function, arising from 
different physical mechanisms, are unrelated. The first contribution (g^ > 0) leads to 
the usual right- and left-hand cuts in the complex energy plane, and does not vanish at 
T = 0. The second piece (g^ < 0) comes from a new cut centered around u = 0, vanishes 
at T = 0, and is unrelated to duality. 

On the hadronic side, at low temperatures T < ^k, only pions from the gas will con- 
tribute to the spectral function. In addition to the (time-like) pion pole contribution to 
/mn|j+\ there is a (space-like) piece in the hadronic spectral function ImU^^^ from the 
center cut in the complex tu-plane. The latter involves an integral of (7r|y4^|27r)(27r|74,^|7r) 

weighted by Bose factors. However, since oc f{ip^ d^i this term appears at order 
T'^/f^ in Imn(s). In addition, the numerical coefficient of this term is very small. Hence, 
a parametrization in terms of the pion pole plus a continuum modelled by perturbative 
QCD should be a good approximation to the hadronic spectral function. 



Using all this information in the FESR Eq.(2.9), suitably extended to T 7^ 0, choosing 
= 0, and neglecting C2 < O2 > (which is proportional to the quark mass squared) one 
obtains the following finite temperature FESR [jl3| 



+ dzM^)[3 - v\z)]nF (^) , (5.4) 

where v{z) = (1 — Am^/z'^)^. Equation (5.4) is an eigenvalue equation fixing so{T) once 
fn{T) is known independently. The advantage of our choice of Green function should be 
evident: at low temperatures (T < fi^) apart from fjriT) there are no other unknown 
T-dependent hadronic parameters such as masses and widths. For instance, QCD sum 
rules for the correlator of two vector currents at T 7^ HiUll) involve Mp{T) and Tp{T) 
which, unlike fniT), are a priori unknown and model-dependent. I shall come back to 
this point in Section 7. 



In the chiral limit the FESR Eq.(5.4) takes the simple form 

Sn'fAT) - f^T^ = f '""^^ ds tank , (5.5) 
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where so{T = 0) = 87rV^(0). In [g the FESR Eqs.(5.4)-(5.5) were solved using the 
low temperature expansion for fTr(T) from 0], i.e. Eq.(3.9). The result is that So(T) 
vanishes at a critical temperature T^i, but < Tc {T^. being the critical temperature 
for chiral-symmetry restoration). However, one should keep in mind that Eq.(3.9) is 
not valid in the vicinity of T^.. As shown in if one uses the expression for fn{T) 
from the composite operator formalism, so(T) vanishes at practically the same temper- 
ature as /7r(T), i.e. — Tc. This is illustrated in Fig. 4 (reproduced from [|]), where 
(so(T)/so(0))^/^ ~ /vr(T)//^(0), except very close to the critical temperature. Given the 
uncertainties of the method, this minor difference can be safely ignored. 



Independent confirmation of this result may be obtained by using the first Weinberg sum 
rule at T 7^ 0, which in the chiral limit is given by 



1 r°° 

- / ds\lmUv{s,T) -lmUA{S,T)] 

TT JO 



(5.6) 



This sum rule was studied in |jT3| using Eq.(3.9), the result being essentially the same as 
with the FESR Eq.(5.5). However, the authors of obtained T^/Tc — 0.99 using the 
more accurate expression for fniT) from the composite operator formalism. 



In summary. Finite Energy QCD sum rules at T 7^ lead to the prediction that so{T) 
vanishes at some critical temperature, which is essentially the same as that for chiral- 
symmetry restoration : ~ T^., provided one uses an expression for f-n{T) valid for all 
T, such as e.g. the one in 0. According to the interpretation of Sq{T) as a relevant order 
parameter for quark deconfinement (see Section 4), one can conclude that the QCD-FESR 
provide evidence for the existence of this phase transition. 



6 Pion form factor at finite temperature 



Independent phenomenological evidence for the deconfinement phase transition in QCD 
may be obtained e.g. by studying the thermal behaviour of the electromagnetic form 
factor of the pion, F^. In this case one would expect the size of the pion to increase with 
increasing temperature. At the critical temperature the pion radius should presumably 
diverge, indicating quark-gluon deconfinement. In this Section I discuss a recent deter- 
mination [20 1 of the T- dependence of F^^ in the space-like region using a Finite Energy 
QCD Sum Rule (FESR). The pion form factor at T = has been extensively studied in 
the past with FESR, as well as with Laplace transform QCD sum rules In order to 
establish some notation, as well as the T = normalization, I briefly describe the method 
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at T = before introducing thermal corrections. 



The appropriate object to study is the three-point function 

n^.A(p,p', <l) = ^'J d'x dS e^(p'^-^^) < 0|T(4(x) VM ^^(0))|0 > , (6.1) 

where A^{x) =: 'u(x)7^75(i(x) : is the axial-vector current, V\ is the electromagnetic 
current, and q = p' — p the momentum transfer. On general analyticity grounds, the 
three-point function Eq.(6.1) satisfies the double dispersion relation 

^,.A(/,p'^Q^) = r ds r ds' (6.2) 

■n^ Jo Jo [s + p^){s' + p'^) 

defined up to subtractions, which are disposed of by Laplace improving the Hilbert trans- 
form, or by considering FESR. The correlator Eq.(6.1) involves quite a few structure func- 
tions, associated with all the Lorentz structures that can be formed with the available 
four-momenta. In principle, it should not matter which particular structure one chooses 



to project the pion form factor. Following pT[ in choosing the combination P^PyP\^ where 



P = p + p', the hadronic spectral function in the chiral-limit reads 

Pis, s', Q')\had = I fMQ')S{s)Sis') + Pis, s\ Q') \qcd[1 - ^(^o - s - s')] , (6.3) 
where sq signals the onset of the continuum, ~ 93 MeV, and 



p{s, s', Q')\qcd = t^^ f3A(x + Q'){x + - - 5Q\x + Q'f] , (6.4) 



to one- loop order (and in the chiral-limit), with 

\ = y^ + Q^i2x + Q^), (6.5) 

and x = s + s',y = s — s'. Since one is interested in writing the lowest moment FESR for 
F^, i.e. 

FUQ^) = ^ r dx r dyp{x,y,Q^)\QCD, (6.6) 

^0 J-x 

rather than a Laplace transform QCD sum rule, the non-perturbative power corrections 
entering the OPE are of no concern here (they contribute to higher moment FESR). The 
integration region in Eq.(6.6) has been chosen to be a triangle in the (s,s') plane, with 
base and height equal to Sq. Other choices of the integration region, e.g. a square region 
of side Si ~ So/^/2, give similar results. The solution to the FESR Eq.(6.6) is 



14 



Although not evident from Eq.(6.7), it is important to reahze that this analysis is only 
valid in the region > 1 GeV^, where one expects a reasonable convergence of the OPE. 
This limitation is of no relevance if one is only interested in the thermal behaviour of the 
ratio F^{Q'^ ,T)/ Ft^{Q'^ ,0). In any case, as shown in Eq.(6.7) provides a reasonable 
fit to the experimental data in the region ~ 1 — 4 GeV^, if Sq — 1 GeV^. 



The spectral function Eq.(6.4) at finite temperature was calculated in |gO| using the 
Dolan-Jackiw formalism. After substitution in Eq.(6.6), the result can be expressed as 

^AQ',T) = —— / dx / dypix,y,Q^)\QCDFix,y,Q^,T), (6.8) 
/2(T) Jo J-x 

with 

F{x, y, Q^, T) = 1 - ni - n2 - ris + nin2 + ^1^3 + n2n3 , (6.9) 



ni=n2 = nF\\^J^^^-Y^\ \ , (6.10) 



_ , , Q^ + (x-y)/2 , , 

^3 = «F F= , (6.11) 

and ni? is the Fermi thermal factor. In the equations above a frame was chosen such that 
= {u, 0), and p'^ = {u', p'), in which case 



- = ,/^ , = (6.12) 

V 2 • 2^ 

It has been explicitly checked that the ratio 

is essentially insensitive to other choices of frames. For instance, one may choose = 
(a;, p), and p'^ = {u', — p), which leads to different arguments in the thermal factors, but 
roughly the same ratio R{T). The temperature dependence of the continuum threshold, 
so(T), is given by (see discussion in Section 5) 



The ratio Eq.(6.13) is shown in Fig. 5 for = 1 GeV^. This result for R{T) is in nice 
agreement with the expectation that as the temperature increases, F^^ should decrease 
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and eventually vanish at the critical temperature for deconfinement T^. 



Although the OPE breaks down at small values of Q^, one may still extrapolate the 
ratio Eq.(6.13) into this region just to study the qualitative temperature behaviour of the 
electromagnetic radius ratio < >t / < >o- Doing this, one finds that this ratio 
increases monotonically with T, doubling at T/Td ~ 0.8, and diverging at the critical 
temperature. This divergence of < >t may be interpreted gnal for quark 

deconfinement. In fact, the behaviour of < >t can be traced back to the temperature 
behaviour of the asymptotic freedom threshold Sq{T). As Sq{T) decreases with increasing 
T, a signature of quark deconfinement, the root-mean-square radius of the pion increases. 
This result is in qualitative agreement with the one obtained in the framework of the 
Nambu-Jona Lasinio model 



7 Laplace sum rules: do they fail ? 



There are a few papers ||2^-|2^ devoted to Laplace transform QCD sum rules at finite 



temperature. I shall concentrate mostly in the vector-vector correlator, e.g. the p^-meson 
channel. All of these analyses suffer from the following serious drawbacks. 

(i) A single Laplace sum rule involves three unknowns: the 
mass and the width of the p- meson, Mp{T) and Vp{T), 
and the continuum threshold sq{T). I am assuming that 
the temperature dependence of the vacuum condensates 
is a known input. No trick of magic allows a determina- 
tion of three unknown quantities from a single equation. 
In principle, at least, one could take derivatives with re- 
spect to the Laplace parameter in order to end up 
with three equations. However, this does not work in 
practice because of (ii) below. 

(ii) The Laplace sum rule is inconsistent with the known 
temperature dependence of the gluon and the quark con- 
densates, as will be shown below. This inconsistency 
translates into a breakdown of the FESR beyond the 
lowest moment. 
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(iii) A zero-width approximation has been used. While it is 
true that at very low temperatures this approximation 
does not differ much from a finite-width parametriza- 
tion, at intermediate temperatures this is not the case. 
Since on physical grounds we expect the width to in- 
crease with increasing T (see Section 4), the zero- width 
approximation is in principle self-contradictory. 



Since the Laplace sum rule is equivalent to an infinite number of FESR, it is simpler to 
break it up into a series of FESR in order to show (ii) above. However, this is not strictly 
necessary; the same conclusion follows from the Laplace sum rule itself, but the procedure 
is a bit more involved. The first three FESR of the type Eq.(2.9) at finite temperature are 
(notice the change of normalization, as I am now discussing the neutral p-meson channel) 

/ ds - ImIio{s,T)\RES = ^ dstanh{—)+—, (7.1) 

Jo 71 Stt^ jo 4T 18 

Ci«Oi»= ds s tanh{^) - 871^ ds s - ImUo{s,T)\iiES , (7.2) 

Jo AT Jo 71 

rso{T) psoiT) 1 

Cg « Oq »= - ds tanhii-) + S^r^ / ds s'^ - ImUo{s, T) \res , (7.3) 

Jo 4T Jo TT 

where Cn « >>= < On > (T), the term C2 « O2 » has been neglected (it 
is proportional to m^), and the hadronic resonant spectral function is given by 

-Imllo[s,I)\RES-^2 (s-M2(T))2 + M2(T)r2(T)- ^^'^^ 



We now have three equations to determine the three unknowns so(T), Mp{T), and Tp{T), 
provided the temperature dependence of the condensates is used as an input. The gluon 
condensate in the chiral limit, and at low T, is given by 

1215/4^ T 4^' ^ ' 

where Ap = 275 ± 65MeV. The above result is valid up to T ~ 100 - 120MeV^. Be- 
yond these temperatures the contribution of massive states, as well as chiral symme- 
try breaking corrections, become gradually more important. Numerically, though, with 
< asG^ >~ O.lGeV^ [|] these corrections are below the 10% level at T ~ ISOMeV 
H. In view of this I shall take the gluon condensate to be temperature- independent. 



Concerning the dimension d = 6 four-quark condensate, I assume it is proportional to 
<< qq and use << qq >> from the analysis of p. At dimension d = 4 and d = 6 

the heat bath can sustain non-scalar condensates Numerically, though, they play a 
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negligible role in this analysis (I have used the estimates of [p7|). 



In trying to solve Eqs.(7.1)-(7.3) I found no meaningful global solution for all three un- 
knowns: So(T), Mp(T), TplT). As I show next, the reason for the absence of a solution to 
all three FESR is that they are inconsistent with the assumed temperature dependence 
of the condensates. 



Let us turn the argument around, and assume a temperature functional form for Mp{T) 
and TplT), and use the three FESR to determine Sq{T), C4 << O4 >>, and Cq « 
Oq ». For definiteness I shall assume that Tp{T) is given by Eq.(4.2). Variations of 
this functional form do not change the conclusions. For Mp{T) I consider two possibil- 
ities: (a): Mp(T) = Mp(0), and (b): Mp{T) as determined using unitarity of the vr — vr 
scattering amplitude, and to lowest order in the virial expansion [0. The latter analysis 
leads to a monotonically increasing Mp{T). Case (a) has already been discussed in p9| ; 
the solutions of the first two FESR for sq{T) and C4 << O4 >>, respectively, are shown 
in Fig. 6. Although sq{T) does decrease with T as expected, the thermal dependence of 
the gluon condensate is definitely wrong. Also, solving the third FESR Eq.(7.3) gives a 
quark condensate that decreases with T even faster than the gluon condensate, again in 
contradiction with expectations. Next, case (b) is qualitatively similar to case (a), as may 
be appreciated from Fig. 7. In this figure, the function Mp{T) is taken from Hence, 
the conclusion from this analysis is that while the first FESR, Eq.(7.1), does lead to a 
reasonable result for sq{T), the higher moment FESR are in contradiction with the known 
T-dependence of the vacuum condensates. Therefore, the FESR program breaks down 
beyond the lowest moment, and so does the Laplace transform. As mentioned before, 
this analysis can be carried out directly with Laplace sum rules, without reference to the 
FESR, but the conclusions are the same. 

There is a third possibility, which I find quite attractive (although it does not solve the 
above inconsistency problem). This is to consider the ratio So(T)/so(0) as a universal 
function, i.e. the same for any channel. Clearly, the value of So(0) will depend on the 
particular channel under consideration, as it is obviously different in the vector and axial- 
vector channels, and the nucleon channel, etc.. However, the ratio so(T)/so(0) may be 
thought of as a phenomenological order parameter for deconfinement, in which case it is 
reasonable to assume it is channel-independent. In this case, one can use e.g. the first 
FESR Eq.(7.1) to predict Mp{T), and the second FESR Eq.(7.2) to predict C4 « O4 ». 
The results of this excercise are shown in Fig.8, where I have used So(T)/so(0) from (see 
Fig. 4). Interestingly enough, the temperature dependence of the p-meson mass turns out 
to be in good agreement with the determination using unitarity of the vr — vr scattering to 
lowest order in the virial expansion [|1^. This supports the validity of the lowest moment 
FESR Eq.(7.1), but unfortunately it does not solve the problem with the higher moment 
FESR or with the Laplace transform, as C4 << O4 >> comes out wrong. Whether 
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additional mechanisms could be brought into play in order to rescue the QCD sum rule 
program remains an open problem. 



Concerning the existing analyses of the vector channel based on Laplace sum rules 
2^ , [p6| , p7[| , they all predict a decrease of the p-meson mass with increasing T. This 
is in contradiction with the result of |19|, which is quite general. In the light of the 
above discussion, it is possible to understand why these analyses predict the wrong T- 
dependence of Mp. Laplace sum rules are inconsistent with the known thermal behaviour 
of the quark and gluon condensates. There is also a Laplace sum rule determination of 
the nucleon mass at finite temperature which gives again a monotonically decreasing 
Mn{T). This is in contradiction with a lowest moment FESR analysis ||T^ which predicts 
an almost constant nucleon mass. It is also inconsistent with [|l4l . The reason behind the 



discrepancy is the same as in the vector channel. I should add that in |^ it is claimed 
that the nucleon mass should decrease, since it is proportional to the quark condensate 
<< qq ». This argument is falacious because the contribution from the hadronic cut 
in the complex energy plane centered around the origin (the so called scattering term) 
prevents the nucleon mass from decreasing [O. 
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